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The nucleation of a droplet of stable cylinder phase from a metastable lamellar phase is examined 
within the single-mode approximation to the Brazovskii model for diblock copolymer melts. By 
employing a variational ansatz for the droplet interfacial profile, an analytic expression for the 
interfacial free-energy of an interface of arbitrary orientation between cylinders and lamellae is found. 
The interfacial free-energy is anisotropic, and is lower when the cylinder axis is perpendicular to the 
interface than when the cylinders lie along the interface. Consequently, the droplet shape computed 
via the Wulff construction is lens-like, being flattened along the axis of the cylinders. The size of 
the critical droplet and the nucleation barrier are determined within classical nucleation theory. 
Near the lamellar/cylinder phase boundary, where classical nucleation theory is applicable, critical 
droplets of size 30-400 cylinders across with aspect ratios of 4-10 and nucleation barriers of 30-40 
ksT are typically found. The general trend is to larger critical droplets, higher aspect ratios and 
smaller nucleation barriers as the mean-field critical point is approached. 
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I. INTRODUCTION 

The decay of a metastable state via the thermally- 
activated formation and subsequent growth of droplets 
of the equilibrium phase is known as nucleation. If the 
nucleation rate is small the formation of critical droplets 
— those droplets large enough to overcome the nucle- 
ation barrier and grow — can be considered to be a 
quasi-equilibrium process. In the classical homogeneous 
nucleation scenario, the critical droplet size and nucle- 
ation barrier are determined from a balance between the 
droplet interfacial and bulk free-energies 0. 

Nucleation from a disordered initial state to a dis- 
ordered final state, as occurs in the liquid-gas transi- 
tion, has been well-studied. The study of nucleation 
when ordered, periodic phases are present is challeng- 
ing, and has been mainly studied in the materials science 
context of crystal growth 0. However there is also a 
large class of soft-materials where competing interactions 
select a microscopic length-scale and lead to ordered, 
periodic microphases Q- Examples include ferrofluids, 
thin-film magnetic garnets, Rayleigh-Benard convection, 
Langmuir monolayers and block copolymers. Whether 
the material is hard or soft, the challenge to modelling 
nucleation in systems with microstructure arises mainly 
from the need to incorporate the various length-scales 
into the theory. In the case of block copolymer melts 
these length-scales are the interfacial width between mi- 
crodomains, the period of the microdomains, the width 



of the droplet interface, and the critical droplet size. 
Another complication is that the symmetry of the un- 
derlying microstructure will lead to anisotropic droplet 
interfacial free-energies and anisotropic critical droplet 
shapes. Furthermore, the droplet interfacial structure 
becomes complicated as the different microdomain sym- 
metries merge into each other. 

In this paper we will focus on diblock copolymer melts, 
whose equilibrium phase behaviour is well-understood 
^ IE 1^- At high temperatures the system is disor- 
dered, while at lower temperatures several ordered peri- 
odic phases — lamellar, hexagonally-packed cylindrical, 
body-centred-cubic spherical, and gyroid — compete for 
stability. This understanding of the equilibrium proper- 
ties of diblock copolymer melts provides a solid founda- 
tion on which to study the kinetics. Within the diblock 
copolymer phase diagram we will study nucleation near 
the first-order lamellar/cylinder phase boundary follow- 
ing a temperature jump into the cylinder phase. This 
is the simplest order-order transition to examine, yet it 
contains all of the challenges and issues mentioned above. 
The goal of this work is to compute the size and shape 
of the critical droplet and the free-energy barrier to nu- 
cleation for a stable cylinder phase nucleating from a 
metastable lamellar background in a diblock copolymer 
melt. 

Previous theories for the kinetics of order-order transi- 
tions in block copolymers have examined the behaviour 
near the order-order spinodal where, in contrast to nu- 
cleation, a particular mode of fluctuation out of the 
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metastable state becomes unstable, and the metastable 
state transforms uniformly throughout the entire sample 
into a more stable state 0, II H 13 IH 111 13 ■ The re- 
cent successful application of self-consistent field-theory 
to block copolymers relies on efficient reciprocal-space 
methods for treating ordered, uniform periodic phases 
IMIMI^' In nucleation, this uniformity is absent, and this 
has limited the application of reciprocal-space methods 
to the problem. The kinetics of various order-disorder 
and order-order transitions have been studied numeri- 
cally in Refs. CI El El E El . The relatively smaU 
three-dimensional system sizes currently accessible nu- 
merically make it impossible to handle the various length 
scales involved in nucleation. Thus direct numerical stud- 
ies are generally limited to examining transitions across 
the entire system at once {i.e. the spinodal regime). It is 
not clear what role, if any, the kinetic pathways and uni- 
form, periodic intermediate states found in these studies 
play in the nucleation of compact droplets of one ordered 
phase in another. 

The existing theoretical work focusing on nucleation 
in polymeric systems considers different systems and 
situations than the one considered here. Wood and 
Wang recently examined nucleation in polymer blends 
using self-consistent field-theory Ell- Fredrickson and 
Binder |2l| and later Hohenberg and Swift 12^ exam- 
ined the nucleation of a lamellar phase from an initially 
disordered phase in melts of symmetric diblock copoly- 
mers. While Ref. only considered a spherical droplet, 
Ref. lU found that the critical droplet shape should be 
anisotropic, with the lamellar normal along the long-axis 
of the droplet. Two-dimensional numerical simulations 
by Nonomura and Ohta have recently examined the dy- 
namics of droplets of lamellae nucleating in a metastable 
cylinder phase . This is close to the situation consid- 
ered here, and we will discuss it in more detail later. 

Early experimental work by Hajduk et al. observed 
a reversible thermotropic transition between the lamel- 
lar and cylinder phase in a polystyrene-poly (ethenc-co- 
butene) diblock copolymer [23. The transition path- 
way from lamellae to cylinders indicated by their data 
suggests that they were observing spinodal decomposi- 
tion rather than nucleation. Sakurai et al. examined the 
reverse transition in poly(styrene-6/ocfc-butadiene-&/ocfc- 
styrene) triblock copolymers, where a non-equilibrium 
cylinder phase produced by a selective solvent transforms 
into a stable lamellar phase upon annealing It ap- 
pears that they were also observing spinodal decompo- 
sition, although the existence of isolated regions of the 
metastable cylinder phase in their system is interesting. 
Floudas et al. have studied kinetics of the lamellar-to- 
cylinder transition in poly(isoprene-Wocfc-ethylene oxide) 
where the ethylene oxide block is crystalline in the lamel- 
lar phase 13 . This crystallinity, the observed rapidity of 
the transition and the lack of specific information about 
the nuclei limit comparison with the present work. Ex- 
perimental work that directly addresses nucleation is lim- 
ited to studies of order-disorder transitions between ei- 



ther the disordered and the lamellar phases |2g or, as will 
be discussed in more detail later, the disordered and the 
cylinder phases IH |2l |3 • These experiments suggest 
that there is an incubation time, followed by nucleation 
of anisotropic droplets of the ordered phase from the dis- 
ordered phase. 

To study the decay of a metastable lamellar state via 
the nucleation of compact droplets of cylinder phase we 
employ the Landau-Brazovskii model H Isl |^. This 
model is appropriate to systems with a weak modulation 
of the order-parameter, or weak segregation, as occurs 
for block copolymers near the mean-field critical point. 
Other phenomena studied using this model include weak 
crystallization, the nematic to smectic-C transition in 
liquid crystals "sT] and Rayleigh-Benard convection 3^. 
The Landau-Brazovskii free-energy can be derived from 
the many-chain Edwards Hamiltonian for diblock copoly- 
mers in the weak segregation limit 0,I3 ■ This important 
connection allows us to move beyond phenomenology in 
this study and make specific predictions about the size 
and shape of the critical droplets, and the nucleation bar- 
rier. 

We apply the single-mode approximation, accurate at 
weak segregation, to the Landau-Brazovskii model, which 
results in an a mplitude m odel studied previously in dif- 
ferent contexts 22,'34','35',''3i^. Since we will work within 
the framework of classical homogeneous nucleation the- 
ory, a key ingredient in our theory is the droplet interfa- 
cial free-energy. We compute the interfacial free-energy 
for a planar interface separating coexisting lamellar and 
cylinder phases from the amplitude model. By using 
a variational approach, we obtain an analytical expres- 
sion for the interfacial free-energy for interfaces of ar- 
bitrary orientation. Previous studies of interfacial free- 
energy, which were primarily numerical, have examined 
only selected interfacial geometries in two dimensions 
[37| . With our results for the interfacial free-energy, we 
can compute the droplet shape from the Wulff construc- 
tion 13 ■ The Wulff construction has been used to study 
the anisotropic droplet shape arising during the disorder- 
to-lamellar transition in Ref. ^21] and more recently in 
Ref. 13- We then apply our results for the interfacial 
free-energies and droplet shape to calculate the critical 
droplet size and nucleation barrier from classical homo- 
geneous nucleation theory. 



II. THEORY 

A. Landau-Brazovskii Theory 

We consider an incompressible melt of n AB diblock 
copolymers in a volume Vq at a temperature T. The total 
degree of polymerization of the diblock copolymer is N. 
The monomer density is thus po = nN /Vq. The degree of 
polymerization of the A block is JaN^ where < /a < 
1. We employ the Landau-Brazovskii theory for weak 
crystallization to study nucleation in this system H IboI 
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l3l| . The position-dependent order-parameter, </)(ro), in 
this theory is defined as the deviation of the normahzed A 
monomer concentration, 0^, from the uniform state, and 
is given by 0(ro) = 4>a{yq) — fA- The Landau-Brazovskii 
free-energy, Fq, is an expansion in terms of this order- 
parameter, 
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In Eq. H2.1|l /o is the Landau-Brazovskii free-energy per 
block copolymer in units of kgT, is the bare correlation 
length, go is the critical wavevector, tq is the reduced 
temperature, and 70 and Aq are expansion coefficients. 
For stability, Ao > 0. The Landau-Brazovskii free-energy 
is able to account for the observed diblock copolymer 
microstructures 9] . 

The Landau-Brazovskii free-energy can be derived 
from the many-chain Edwards Hamiltonian for diblock 
copolymers following the method of Leibler j^] and Ohta 
and Kawasaki ;33j. Like many authors, we have assumed 
that the third- and fourth-order expansion coefficients 
are local. The truncation of the order-parameter expan- 
sion at fourth-order in Eq. H2.1|l and the assumption of a 
single dominant wavevector, q^, in the gradient term are 
valid if the system is in the weak-segregation limit near 
the disorder-to-order spinodal and the mean-field critical 
point. Since, in mean-field theory, the lamellar-cylinder 
transition occurs in this region, the Landau-Brazovskii 
free-energy is appropriate to our study. 

It is convenient to rescale Eq. (|2.1|) by expressing 
lengths in units of Qq^ and the free-energy is in units 
of Ao. Under the rescalings 
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the Landau-Brazovskii free-energy becomes 
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One of the goals of this work is to calculate the critical 
size, shape, and nucleation barrier for droplets of sta- 
ble cylinder phase nucleating from a metastable lamellar 
phase in terms of experimentally measurable quantities. 
To do this we must relate the parameters qq, ^q, tq, 70, 



and Ao appearing in Eq. H2.1|l to such quantities. By fol- 
lowing the derivation of Eq. (|2.1I) in Ref. |3^ , we can ex- 
press these parameters in terms of x, Ja and Rg that 
appear in the many-chain Edwards model. Here x is 
the Flory-Huggins interaction parameter characterizing 
the repulsion between A and B monomers. The diblock 
copolymer radius of gyration, i?g, is defined through 
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where h is the Kuhn statistical segment length. We have 
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The notation follows that of Leibler In particular, x* 
is the position of the minimum of the function F{x, Ja) 
appearing in the scattering function of Ref. The 
disorder-to-order spinodal, (xA^)s, and the quantity, c, 
are defined through 
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Finally, the vertex functions Fa and F4(0,0), which are 
functions of fA , are computed in Ref. [1| . In Table U 
we reproduce the values tabulated by Fredrickson and 
Helfand for these mean-field parameters in the weak- 
segregation region "ZSj. We will take Ja < 0.5 in this 
paper without loss of generality. 



TABLE I: Mean-field model parameters in the weak- 
segregation region, from Ref. [i^. 



fA 


X* 


{xN)s 


c 


ivr3 


ivr4(o,o) 


0.50 


3.7852 


10.495 


0.4812 


0.0 


156.56 


0.45 


3.7995 


10.698 


0.4844 


-8.608 


169.19 


0.40 


3.8433 


11.344 


0.4945 


-18.81 


212.22 



In terms of these parameters, the scaled quantities ap- 
pearing in Eq. fi.Sfl are 
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For example, for /a = 0.5 we have = 0.0440 and 
7 = 0, while for = 0.45 we have = 0.0413 and 
7 = 0.0509. When xN = 11.0, r = -6.45 x 10^^ ^j. 
/a = 0.5 and r = -3.56 x lO^^ for Ja = 0.45. The free- 
energies of the equilibrium lamellar and cylinder phases 
computed in the weak-segregation regime using these pa- 
rameter values differ by only a few percent from state-of- 
the-art self-consistent field-theory calculations |4l|. 



B. Geometry of the Droplet and the Single Mode 
Approximation 

In principle, to describe a critical droplet of stable 
cylindrical phase in a background of metastable lamel- 
lar phase one needs to find a saddle-point solution to the 
Euler-Lagrange equations corresponding to Eq. (|2.8() . In 
practice, for the three-dimensional, anisotropic droplets 
we consider, finding such a solution by a direct numerical 
method is not, at present, feasible due to the large system 
sizes required to resolve the various length-scales in the 
problem — the A/B interfacial width, the microstruc- 
ture period, the droplet interfacial width and the droplet 
size. Thus, to make progress, we need to make some ba- 
sic assumptions about the geometry of the droplet and 
the form of the order-parameter. 

We consider a droplet of stable cylindrical phase in a 
background of metastable lamellar phase with the geom- 
etry shown in Fig. ^ Inside the droplet, the cylinders 
are ordered into a hexagonal lattice and aligned along 
the z-axis. Outside the droplet, the lamellar normal is 
directed along the y-axis. We assume an epitaxial rela- 
tion between the cylinders and the lamellae. Matsen's 
work on the cylinder/gyroid and cylinder/sphere tran- 
sitions shows that these order-order transitions proceed 
approximately epitaxially, and that imposing an epitaxial 
relation results in only a minor change in the energetics 
of the transition 0, 0| . This suggests that the epitax- 
ial assumption may also have only a slight effect on the 
energetics of the lamellar/cylinder transition. Further- 
more, experiments indicate that at the lamellar/cylinder 
transition the cylinders will form with their axis in the 
plane of the layers, as we have assumed (23|, i2j|. The 
lamellar period, Di, is related to qq by 

27r 

go = j^. (2.20) 

The distance between cylinders, is epitaxially related 
to Di by Dc = 2Di/V3. 

The orientational epitaxy enables us to describe both 
the cylindrical and the lamellar order in terms of a single 
set of reciprocal-lattice vectors , indexed by the integer 
i. Since we work in the weak-segregation regime it is 
sufficient to restrict ourselves to the first mode, those 
vectors with |Gi| = qo (in scaled units, |Gi| = 1), instead 
of using the complete set of vectors G.;. In the single- 
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FIG. 1: Schematic picture of the droplet geometry. In reality 
the critical droplet is much larger than shown, (a) View of 
the x~y plane. The lamellar normal is along the y-axis and 
the cylinders form a hexagonal lattice inside the droplet with 
their axes along z (out of the page). The lamellar period is Di, 
the distance between cylinders is D^. The azimuthal angle, 
0, is indicated, (b) Perpendicular view to (a), showing the 
orientation of the cylinders relative to the droplet interface. 
The polar angle, 9, is indicated. 

mode approximation, the order-parameter is written as 

0(r) = 2ai(r) cos(Grr)+2a2(r) [cos(G2-r)+cos(G3-r)], 

(2.21) 

where oi and 02 are spatially-dependent amplitude func- 
tions which define the droplet geometry and the Gi are 
given (in scaled units) by 

Gi = y (2.22) 

G2 = i(-V3x + y) (2.23) 

G3 = i(-V3x-y). (2.24) 

Pure lamellar order is described by ai 7^ and a2 — 0, 
pure cylindrical order by oi = 02 7^ 0. 

If we assume the amplitudes ai and 02 are slowly vary- 
ing on the scale of Dc, we can separate the length-scale for 
variations in the amplitude, the droplet interfacial width, 
from the length-scale of the underlying microstructure. 
Within this slowly- varying amplitude approximation, the 
Landau-Brazovskii free-energy, Eq. (|2.8|) , can be written 



in terms of the amplitudes as 

+ (G2 • Vaa)' + (G3 • Vaa)'] + r(a? + 2al) 
—270102 + ^(fli + + Sajflj) |. 
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Extremization of this free-energy produces two Euler- 
Lagrangc equations for the amphtudes, 

2S_^V^ai~8C^{Gi-Vfai+vi{ai,a2) = (2.26) 
4e2vV-8e'[(G2- V)2 + (G3- V)2]a2 + «2(ai,a2) 

= 0, (2.27) 



where 



wi (01,02) — 2rai — 27O2 + ( 
^2(01, 02) = 4ra2 — 4701O2 + 603 + 4oi02. (2.28) 



The free-energy H2.25|l and the amphtude equations 
(|2.26|) and H2.27|l form the basis for our analysis of the 
droplet nucleation problem. 



C. Equilibrium phase diagram 

Uniform, periodic lamellar and cylindrical states are 
position- independent solutions of Eqs. (|2.26l) and H2.27|l . 

When oi = 02 = 0, the system is disordered. The uni- 
form lamellar solution, which exists for t < 0, is 
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(2.29) 
(2.30) 



From Eq. (|2.25(l we see the uniform lamellar solution has 
free-energy 



(2.31) 



The uniform cylindrical solution has oi = 02 = Oc with 
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and has a free-energy 



fc = 3Ta~ - 270; 



(2.32) 



(2.33) 



For 7 > the solution with the positive root in Eq. (|2.32() 
has the lower free energy and corresponds to cylinders of 
A in a B matrix. The case for 7 < can be constructed 
from the 7 > case by recognizing that the free-energy, 
Eq. (|2.33() . is invariant under 7 — > — 7 and Uc —> —Uc- It 
is straightforward to show that, along a phase boundary 
or spinodal, the relations r — xj'^ and Oc = T^c hold for 
7 of either sign, with 



1 + VI - lOx 



(2.34) 




FIG. 2; Equilibrium phase diagram for the Landau- 
Brazovskii model in the single-mode approximation, Eq. 
12.2511 . The solid lines are phase boundaries between the 
indicated phases. The dashed line is the disorder-to-order 
spinodal. The dotted line is the lamellar-to-cylinder spin- 
odal. The dashed-dotted line is the cylinder-to-lamellar spin- 
odal. The body-centred-cubic spherical phase, which is stable 
between the cylinder and disordered phases, and the gyroid 
phase, which self-consistent field-theory predicts is stable for 
XN > 11.14 along the lamellar/cylinder phase boundary, are 
not shown. The second-order, mean-field critical point is in- 
dicated by the dot. 



and X a constant. The lamellar/cylinder phase boundary 
occurs for 



7 + 3V6 



^2.8697, 



(2.35) 



as noted in [33|, and the cylinder/disorder phase bound- 
ary occurs for x = 4/45. The stability limits, or spin- 
odals, for these uniform phases can be found by expand- 
ing Eq. I|2.25|l to second-order in deviations from the uni- 
form solution and looking for parameters for which the 
matrix of partial second-derivatives has a negative eigen- 
value. The disorder-to-order spinodal is just r = 0. On 
increasing temperature (increasing r, decreasing x^) 
stability limit for the lamellar phase is reached when 



-27^ 



(2.36) 



At this point the lamellar structure will spontaneously 
transform into the cylinder phase. The spinodal for the 
reverse, cylinder-to-lamellar, transition occurs when r = 
-872. 

We construct the equilibrium phase diagram for the 
free-energy Eq. (|2.25|) in Fig. [3 in terms of x^ and /a, 
using relations H2.18|l and (|2.19|l . In mean- field theory 
there is a second-order critical point at /a = 0.5 and 
xN = 10.495 (7 = r = 0), as indicated. In self-consistent 
field-theory calculations the body-centred-cubic spheri- 
cal phase is stable between the cylindrical and disor- 
dered phases and the gyroid phase is stable along the 
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lamellar/cylinder phase boundary for > 11-14 and 
/a < 0.452 {fA > 0.548) 5]. We have not included these 
phases in our diagram since our primary concern is with 
the lamehar/cyhnder transition and since we work near 
the lamellar/cyhnder phase boundary for /a > 0.45. As 
mentioned previously, the free-energies of the equilibrium 
phases calculated from Eq. (|2.25|) are within a few per- 
cent of those obtained from self-consistent field-theory. 

D. Planar interface between the lamellar and 
cylindrical phases 

For a planar interface separating the lamellar and 
cylindrical phases the amplitudes ai and 02 become 
functions of the perpendicular distance to the interface, 
s — n • r, only. The unit normal to the interface is 
n = (cos 0sin0, sin (/)sin 0, cos 0), expressed in terms of 
the angles defined in Fig. ^ From Eq. H2.25|l . the free- 
energy per unit interfacial area, /, is then 

/ = J ds{e{a'ir +2e{4)' +miGi-nna[r 

+{(G2 • n)2 + (G3 • 6)^(4)'] + rial + 2al) 
~2-iaxal+'^{a\ + &a2 + Sa\al)^, (2.37) 

where primes denote differentiation with respect to s. We 
note that 

(fi- Gi)^ = sin^^ siT?e (2.38) 
(fi- Ga)^ + (n- G3)2 = i(l + 2 cos^ 0) sin^ 0(2.39) 

Thus terms involving the squares of single derivatives of 
the field (square-gradient terms) do not contribute to / 
when the interface normal is in the z direction {6 — Q). 

III. ENERGETICS OF THE 
LAMELLAR/CYLINDER INTERFACE 

The aim of this section is to compute the interfacial 
free-energy for a planar interface of arbitrary orienta- 
tion n between coexisting lamellar and cylindrical phases. 
Once the interfacial free-energy is known, the Wulff con- 
struction [3^ can be used to find the minimal-energy 
droplet shape, as we show in the next section. It is suffi- 
cient to approximate the interfacial free-energy near the 
lamellar/cylinder phase boundary by its value at bound- 
ary since the interfacial free-energy function is continuous 
through the transition. 

A. An analytic calculation of the interfacial 
free-energy 

Rather than numerically solving the Euler-Lagrange 
equations resulting from Eq. (|2.37|) to get the interfacial 



profile and then computing the interfacial free-energy, 
we instead employ a variational ansatz for the interfa- 
cial profile which allows the interfacial free-energy to 
be calculated analytically. Our variational solution to 
the problem will, of course, produce an interfacial free- 
energy which is greater than that obtained through an 
exact solution of the Euler-Lagrange equations. How- 
ever, our analytical solution will give us much insight 
into the behaviour of the interfacial free-energy along the 
lamellar/cylinder phase boundary. 

We employ the following variational ansatz for the am- 
plitude profile of a planar interface between coexisting 
cylinder and lamellar phases. 




The function h has the properties that h{s) ±1 as s ^ 
±00 and h{0) = 0, thus it describes a pure cylinder phase 
for s — > — 00 separated by an interface at s = from a 
pure lamellar phase, obtained as s — > 00. The form of 
h will be specified later. The variational parameter, w, 
characterizes the interfacial width, and will depend on 
the interface orientation. Substitution of Eqs. H3.1|l and 
(|3.2() into the free-energy, Eq. (|2.37l) . and rescahng s to 
extract factors of w results in the excess free-energy due 
to the interface, 

f - Jh = 5 \ h 7 33 w, 3.3 

where 

h = Jds fi (3.4) 

is the free-energy contribution from the two uniform 
phases (/; = fc at the lamellar /cylinder phase bound- 
ary). The coefficients in Eq. H3.3|l are 



where 



91 = 




- 2aidc + 3al) 


(3.5) 
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+ 2ai{ai — 2ac) sin^ </>] sin^ 9 


(3.6) 
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^ — 00 


x{a1 + 20,2) — 2aia2 






+i(atH 


- 6dj + Safal) + x'^Y 


(3.7) 




h 


= ds {h"f 

J —00 


(3.8) 




h 


nOO 

^ ds [h' f. 

J —00 


(3.9) 



The quantity di = ^/—2x, and and x are given by Eqs. 
(I2.34|) and H2.35|) . In Eq. H3.7|l the functions di and 02 are 
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given by Eqs. (|3.1(l and (|3.2() with ai and Oc replaced by 
ai and respectively, and with u; = 1. The coefficients 
(/I, 52, and are positive. The excess free-energy, Eq. 
is minimized when 



u;2 = (i^*) 



(3.10) 



Substitution of this value for w in Eq. I|3.3|l provides the 
interfacial free-energy, a{9,4>) = /* — //j, for a planar 
cylinder/lamellar interface with normal n, 



i T9\ ^ i T9-2 



4 * 

7 93W 



(3.11) 



The angular dependence in Eq. 1)3.11(1 is contained in 
g2 and w*. Equation H3.6|l indicates that the interfacial 
free-energy depends on the polar angle 9 through factors 
of sin^ 9. Equation (|3.6() also shows that the interfacial 
free-energy is invariant under (f) (j) + nn, where n is an 
integer (two- fold azimuthal symmetry). 

To examine the behaviour of a{9, 0) along the phase 
boundary near the mean-field critical point, where 7 — > 0, 
we have to distinguish between two cases, depending on 
the relative values of 9 and 7. For I7I ^ 32, as obtains 
for finite 9 and 7 ^ 0, we have 



w 



0(|7l) 



(3.12) 



and 



a{9,ct>)=2^h\'y^^+Oi\jf). (3.13) 
For 9-^0 such that 92 ^ I7I we have 

£.92 



w 



l7l 



53 



1 



and 



= |C^/^l7r/^(35i5l)^/' 



2|7l(125i53)i/2 

(3.14) 

3^52 
2|7|(12gi53)i/2 

(3.15) 
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The width of the interface diverges at the mean-field crit- 
ical point, and the interfacial free-energy vanishes there. 
It is interesting to note that the interfacial free-energy 
remains anisotropic as 7 ^ 0. Since 7 ~ -v/r along the 
phase boundary, the scaling with 7 seen in Eqs. ((3.12|l 
and H3.13|l corresponds to the expected mean-field re- 
sults, w-- and a ^ r^/^ [12 . When 32 < I7I 
Eqs. (|3.14() and (|3.15() show that, as 7 ^ 0, the interfa- 
cial width diverges more slowly than the mean-field result 
and the interfacial free-energy vanishes more rapidly than 
the mean-field result. As 7 ^ 0, the condition 52 ^ I7I 



restricts the region over which this behaviour holds to 
very small 9. The non-mean-field scaling in Eqs. (|3.14|) 
and (|3.15|) arises because, for 9 = 0, only the square- 
second-derivative terms in Eq. (|2.37l) are non- vanishing. 
Finally, far from the critical point, as I7I — > 00, the con- 
dition 52 ^ 7 holds for all 9 and the expressions (|3.14|) 
and (|3.15|) apply. At leading order w* ^ and cr ^ 00 as 
I7I — > 00. The interfacial free-energy becomes isotropic 
in this limit. 

None of these results depend on the form chosen for h 
in the variational approximation. The only approxima- 
tion in Eqs. 1(3. l|l and 1(3. 2|l is that the interfacial profiles 
for two different orientations have the same basic shape, 
h, and can be scaled onto one another by an orientation- 
dependent scale-factor w. 



B. Results for the lamellar /cylinder interfacial 
free-energy 

To calculate the interfacial free-energy we will now 
choose h{s) = tanhs, which is a reasonable approxima- 
tion for an interfacial profile, and satisfies the boundary 
conditions mentioned after Eqs. (|3.1|) and H3.2|l . With 
this choice Ii = 16/15, h = 4/3 and 53 w 0.716 286. 
The unsealed (i.e. actual) interfacial free-energy, do, is 
related to a in Eq. (|3.11l) by 



An 



CTO = 



po b ksT 

7V1/2 



(3.16) 



In Fig. 121 we present results for ctq at the phase bound- 
ary, in units of po b ksT/N^/'^, for Ja = 0.45. Figure 
|31 qualitatively represents the features of the interfacial 
free-energy for all Ja near the mean-field critical point. 
The overall magnitude of the interfacial free-energy in- 
creases with stronger segregation. In Fig. 13^- we show 
CTo as a function of for = 0. When the cylinder 
axis is perpendicular to the interface [9 — 0) the in- 
terfacial free-energy is about 4 times smaller than when 
the cylinders are parallel to the interface [9 = tt/2). In 
Fig. ISJj the interfacial free-energy as a function of (j) in 
the x-y plane [9 = 7r/2) is shown. The interfacial free- 
energy is highest when the lamellae lie perpendicular to 
the interface (0 = 0) and lowest when the lamellae are 
parallel to the interface (</) — 7r/2). The variation of 
the interfacial free-energy with (f) is less than its vari- 
ation with 6*, since 2ai{ai — 2dc) ^ 35,^ in Eq. (|3.6|l . 
Although the interfacial free-energies for — tt/6 and 
4> = 57r/6 are equal, as they should be, it is somewhat 
surprising that the presence of cylinder lattice planes at 
these angles does not produce any features in the inter- 
facial free-energy. It may be that at higher segregation, 
where one has to go beyond the single-mode approxima- 
tion, the existence of cylinder lattice planes will become 
manifest in the droplet shape. Similar work by Netz et 
al. examined the lamellar/cylinder interfacial free-energy 
numerically when 9 = (f> = 7r/2 j37j . They found the 
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FIG. 3: Interfacial free-energy, ctq in units of po b ksT/N^^^ , 
of the lamellar/cylinder interface at the lamellar/cy Under 
phase boundary for /a = 0.45 (x-^ = 11-3263) as per Eqs. 
(ISTTTl and llO^ . (a) As a function of 9, for (/> ^ 0. (b) As 
a function of (j) for 6 — n/2 (in the x~y plane). The angles 9 
and (j) are defined in Fig. 



mean- field scaling, (Tq '--^ r'^/^, near the critical point in 
their model. Self-consistent field-theory has been used 
recently to study the energetics of kink- and twist-grain 
boundaries in diblock copolymer melts — situations also 
involving complicated interfacial structures Q| • 



IV. WULFF CONSTRUCTION OF THE 
CRITICAL DROPLET 

In general, the nucleating droplet of cylindrical phase 
is not spherical and it is necessary to calculate the droplet 
shape from the anisotro pic interfacial free-energy, using 
the WulfF construction Once the droplet shape is 



known the free-energy and size of the critical droplet can 
be calculated. 



A. Wulff construction for the droplet shape 

The droplet shape is found by minimizing the droplet 
surface free-energy, subject to the constraint of constant 
droplet volume. Thus we minimize the function 



Wulff 



drop 



2fiVdrc 



where is a Lagrange multiplier, Vdrop 
droplet volume and 



drop 



dS a{n) 



(4.1) 
is the fixed 

(4.2) 



is the integral of the interfacial free-energy, Eq. H^.llfl . 
over the droplet surface. Implicit in Eq. H4.1|l and in the 
classical nucleation theory we employ is the assumption 
that the droplet interface width is negligible compared 
to the droplet size, so that a separation may be made 
between the bulk and surface of the droplet. We will 
show this to be the case a posteriori, close to the lamel- 
lar/cylinder phase boundary. 

The minimization of Eq. H4.1|) is performed by choos- 
ing, for example, the x-y plane and characterizing the 
shape of the droplet by the height of the interface, 
hoi'^±), above this plane at x_l = {x,y). Minimization 
then leads to the following formula for the height [45l |. 



/i(x_l) = [g(m) m • x_l] 



mm m' 



where 



m 



L 

T 



5(m) - (1 + m^) 



2M/2 g^["(m)] 



(4.3) 

(4.4) 

(4.5) 
(4.6) 
(4.7) 



We have written /i = o'max/L with <7max = 
[(T(fi(m))]max m and L a length-scale determined by the 
value chosen for Vdrop- Equation 14.3|l is the essence of 
the WulfT construction. The minimization in Eq. 14.3(1 
over the variable m for a given x^ is performed numeri- 
cally. 

The droplet shape along the lamellar/cylinder phase 
boundary for /a = 0.45 and 0.49 is shown in Fig. 01 In 
this figure we have chosen the unit of length to be the 
largest dimension of the droplet, in this case the half- 
length, Ix, along the x-axis. The other dimensions are 
determined, in these units, through the Wulff construc- 
tion. Figure shows that the droplet is lens-shaped, 
being flattened along the axis of the cylinders (z-axis). 
This follows from the lower interfacial free-energy for in- 
terfaces perpendicular to the cylinders, compared to in- 
terfaces parallel to the cylinders. The anisotropy of the 
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FIG. 4: Cross-section of the droplet shape calculated via the 
Wulff construction. We have chosen the unit of length to be 
the largest dimension of the droplet, in this case the half- 
length along the a:-axis. (a) Cross-section in the x-z plane for 
/a ~ 0.45 (outermost curve) and /a ~ 0.49 (innermost curve) 
along the lamellar/cylinder phase boundary. The cylinders 
inside the droplet are oriented along the z-axis (vertically) 
and the lamellae outside the droplet lie in the plane of the 
page, (b) Cross-section in the x-y plane for /a = 0.45 (out- 
ermost curve) and /a = 0.49 (innermost curve) along the 
lamellar/cylinder phase boundary. These curves are essen- 
tially indistinguishable. The cylinders inside the droplet are 
in a hexagonal lattice, with their axes pointing out of the 
page. The lamellae outside the droplet are oriented horizon- 
tally, with their normals along the y-axis. 



droplet increases as the mean-field critical point is ap- 
proached. The droplet shape anisotropy can be char- 
acterized by the ratio of Ix to the half-length along the 
z-axis, Iz. From Eqs. H4.3|l and 14.7|l this aspect ratio is 



which, from Eqs. (|3.13() and H3.15|l . scales as 



(4.9) 



as the critical point is approached, and approaches 1 as 
I7I 00. The increase in the aspect ratio of the droplet 
as the mean-field critical point is approached is is a con- 
sequence of the gradient structure of the theory, which 
leads to different scaling of the interfacial free-energy 
with I7I in Eqs. and (|XT5|) . 

Figure 2|3 shows that the droplet is also slightly flat- 
tened along the y-axis, compared to the oi-axis. This 
reflects the trend seen in Fig. where the interfacial 
free-energy is lower when the lamellae lie along the inter- 
face than when they are perpendicular to it. As expected, 
the droplet shape anisotropy is weaker in the x-y plane 
than it is in the x-z plane. The change in the droplet 
shape anisotropy in the x-y plane with Ja is also very 
small. If ly is the half-length of the droplet along the 
y-axis, Eqs. H4.3|l . (|4.7|1 . and H3.13|l lead to 



/ g2(V2,7r/2) 
52(^2,0) 



= 0.9074, 



(4.10) 



^(V2,0) 
fT(0,0) ' 



(4.8) 



at leading order in I7I as the critical point is approached. 
Thus near the mean-field critical point the droplet will 
become extremely flattened along the z-axis, with a 
slightly non-circular cross section in the x-y plane. As 
I7I 00, Eq. (|3.15|) indicates that the droplet shape 
will become isotropic, and ly/lx 1- As anticipated, 
the presence of cylindrical lattice planes at (/) — n/G and 
(j) — 57r/6 is not manifest in the droplet shape in the x- 
y plane. We will compare these results for the droplet 
shape with relevant experimental and numerical results 
in the Discussion section. 



B. The critical droplet 

The droplet volume serves as a reaction coordinate 
during the transformation from the metastable lamellar 
phase to the stable cylinder phase. Once this volume 
is selected, the droplet shape that minimizes the surface 
free-energy is determined from the WulfF construction. 
We will employ the classical homogeneous nucleation the- 
ory P] , in which the critical droplet corresponds to a total 
free-energy maximum along the volume reaction coordi- 
nate, and the total free-energy is a sum of surface and 
bulk terms. The computed droplet shape is independent 
of the value chosen for the droplet volume, thus the shape 
of the critical droplet will be just that calculated in the 
last section. 

When we choose the unit of length to be the largest di- 
mension of the droplet, as in the last section, we can com- 
pute the droplet volume, Vdrop, and surface free-energy, 
Sdrop, in these units as functions of fA- When /a = 0.45, 
for example, Sdrop = 1-949 57 xlO"'^ and Vdrop = 0.784 



000 (the interfacial free-energy in Sdrop, Eq. (|4.2(l . is di- 
mensionless). If we wish to express in units of Qq^, 
the surface free-energy S{lx) and droplet volume V{lx), 
in units of Qq^ and q^^ respectively, will scale as 

S{1.) = Sdrop ll (4.11) 
ViL) = Vdropll (4.12) 

for arbitrary l^- 

We now examine the system slightly away from the 
phase boundary, making the standard assumption of clas- 
sical nucleation theory that the interfacial free-energy 
and interfacial width do not change significantly near the 
phase boundary. When the lamellar phase is metastable 
and the cylinder phase is stable we have 

A/ = < 0. (4.13) 
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Separating the free-energy Fq (in real free-energy units) 
into a bulk and a surface term, and subtracting the free- 
energy, Fi , of the metastable uniform lamellar phase we 
have 



(Fo - Fi)N Ao 



4 



[Vdrop ll ^ + Sdrop ll] ■ (4.14) 



The critical droplet occurs for the half-length, Ic 
maximizes this expression, namely 



Ic 



^Sdro 



SVdrop^f 



for which the nucleation barrier, AF^. = F^rit 
given by 



Ao 



PO 



h^N^^kBT (6x*)3/2 27{VdropY{Aff 



that 



(4.15) 



F. 



(4.16) 



The critical half-length in Eq. H4.15|l is expressed in units 
of (/q"^. Equation (|2.10() has been used in Eq. (|4.16() . 

In Fig.|Slwe plot AFc in units of pob^N^^'^ksT near the 
lamellar/cylinder phase boundary for Ja = 0.45 and 0.49. 
The nucleation barrier diverges at the phase boundary, 
but remains finite at the spinodal — behaviour to be 
expected from classical nucleation theory. Physically, the 
nucleation barrier should go to zero at the spinodal, but 
this is beyond the classical nucleation approach employed 
here. Over the region in from the spinodal to the 
phase boundary the general trend is for the nucleation 
barrier to decrease as the critical point is approached. 
The importance of Fig.|Slis that it allows one to calculate 
the magnitude of the nucleation barrier in units of fc^T, 
once N, po and b are known. Since AFc scales as N^^^, 
we expect polymers with higher indices of polymerization 
will have larger nucleation barriers. Finally, in Fig. [S] we 
plot the nucleation barrier obtained using our theory for 
the interfacial free-energy, but with a spherical droplet 
shape (dashed curve). In this case AFc is significantly 
increased, indicating the importance of using the proper. 



(b) 



Z 1.0- 

O 

§ 0.5- 
c 



<l 




520 10.522 10.524 10.526 10.528 



FIG. 5: Nucleation barrier (solid curves) as a function of xN 
near the lamellar/cylinder phase boundary, (a) For /a ~ 0.45, 
where the phase boundary occurs at = 11.3263 and the 
spinodal occurs at x^ = 11.1358. (b) For /a ~ 0.49, where 
the phase boundary occurs at x^ ~ 10.5285 and the spinodal 
occurs at x^ = 10.5207. For comparison, the dashed curves 
are the result of using the anisotropic interfacial free-energy, 
Eq. lITTTl . but assuming a spherical droplet shape. 



anisotropic droplet shape for an accurate calculation of 
the nucleation barrier. 

In Fig. we examine the dimensions of the critical 
droplet, lex and Icz obtained using Eqs. H4.8|) and H4.15|l . 
as a function of xiV near the lamellar/cylinder phase 
boundary for /a = 0.45 and /a = 0.49. In this fig- 
ure the lengths are in units of the cylinder spacing Dc- 
In addition to lex and Icz we also show the droplet in- 
terfacial widths, Wx and Wz, obtained from the calcula- 
tion of the interfacial free-energy at the phase boundary 
for interfaces oriented with normals along the x- and 
2;-axes, respectively. The droplet size diverges as the 
phase boundary is approached, as expected from classi- 
cal nucleation theory. In Fig. we examine the region 
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FIG. 6: Critical droplet dimensions, in units of the distance 
Dc between cylinders, as a function of near the lamel- 
lar/cylinder phase boundary, (a) /a = 0.45 (b) /a — 0.49. 
The upper solid curve corresponds to lex, the lower solid 
curve to Icz- The interface widths, calculated at the lamel- 
lar/cylinder phase boundary, are indicated: Wx (dashed line), 
Wz (dotted line). The transitions and spinodals occur in the 
same places as in Fig. |3 



in from the spinodal to the phase boundary. The 
general trend is for the scale of the critical droplet in 
this region to increase as the mean-field critical point is 
approached. As anticipated from Eqs. (|3.12|) and (|3.14|) 
the width of the droplet interface increases as the critical 
point is approached. Netz et al. also observed widening 
of the lamellar /cylinder interface as the segregation de- 
creased 113 ■ For these /a and we find lex > Wx and 
Icz > Wz- For /a =0.49 there is a region of xiV for which 
Icz < Wx- In the next section we will discuss the region 
of validity of our theory, and the range of droplet sizes 
and nucleation barriers expected in this region. 



V. DISCUSSION 

To satisfy the assumptions of the Wulff construction 
and the classical nucleation theory we study the region 
close to the lamellar/cylinder phase boundary, but not 
so close that nucleation is rendered unobservable due to 
extremely high barriers. In small-molecule binary fluid 
mixtures, Cahn and Hilliard estimated that a free-energy 
barrier below 60 ksT should produce an observable rate 
of nucleation (at least one nucleation event per day per 
cubic centimetre, say) 0- Since the relaxational dynam- 
ics in polymeric fluids is slower than in small-molecule 
fluids, we expect that a somewhat smaller free-energy 
barrier is required for nucleation to be observable in our 
system. However, in the absence of a kinetic theory for 
nucleation in polymers, we will take 60 ksT to be an 
upper limit for the nucleation barrier. In the context of 
binary polymer blends. Binder and, more recently, 
Wang [43 have argued that mean-field theory breaks 
down when the free-energy barrier becomes less than 10 
ksT- Following these authors, we take a nucleation bar- 
rier of 10 ksT to be the lower limit for applicability of 
the present theory. 

To be consistent with the slowly- varying amplitude ap- 
proximation, the width of the droplet interface should be 
larger than the microstructure period. For Ja — 0.45 
we have Wx ^ Wy Dc, while closer to the mean-field 
critical point Wx ^ Wy > Dc- We typically find Wz < Dc 
(when /a = 0.45, Wz ~ 0.41 D^ and it is not until 
Ja = 0.492 that Wz = Dc)- However, since the vari- 
ations in the underlying microstructure are in the x-y 
plane, relatively rapid variations of the amplitude in the 
z-direction should not affect the validity of the slowly- 
varying amplitude approximation. For /a = 0.45 and 
XN = 11.25 with N = 1000 and po = b^^ the nucle- 
ation barrier is 39 fc^T and the critical droplet, which 
has an aspect ratio of 4.3, is about 30 cylinders across. 
These numbers are reasonable, so we take Ja ~ 0.45 as 
the lower bound on /a for which the theory is valid. For 
Ja < 0.452 self-consistent field-theory predicts that the 
gyroid phase is stable along the lamellar/cylinder phase 
boundary. Furthermore, the weak-segregation approx- 
imation will become increasingly inaccurate as /a de- 
creases. 

As the mean-field critical point is approached the trend 
is for the nucleation barrier to decrease, the droplet in- 
terface width to increase, and the droplet size to in- 
crease. Thus the slowly-varying amplitude approxima- 
tion should be increasingly accurate as the critical point 
is approached. We also have lex > Wx and Icz > Wz in 
this limit, which is a necessary condition for the classi- 
cal nucleation and Wulff approaches to be valid. It is 
possible for Icz < Wx for some super-heatings, as seen 
in Fig. IHb, possibly indicating a breakdown of this ap- 
proach, however one can go to smaller super-heatings 
where Icz > Wx- As an example, for Ja = 0.49 and 
xN = 10.5275 with N = 1000 and po = 6"^ the nucle- 
ation barrier is 28 fc^T and the critical droplet, which 
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has an aspect ratio of about 10, is about 420 cylinders 
across. For these values of /a and find l^z ~ 3.7 

Wx- The interfacial width for interfaces parallel to the 
cylinders is on the order of 5-6 D^, while for perpendic- 
ular interfaces the width is on the order of 0.9 Dc- Near 
the mean-field critical point fluctuations will be impor- 
tant and will renormalize the basic model, Eq. H2.1|l . as 
discussed in 21, 30, 40). Mean- field theory is recovered in 
the limit iV — > oo. The technique discussed here can be 
applied to study nucleation in the renormalized model, in 
which the lamellar/cylinder phase boundary terminates 
at a lamellar/cylinder/disorder triple point for Ja < 0.5, 
where the nucleation barrier is expected to remain finite 
|40| |. In addition, we expect that when the nucleation 
barrier becomes less than about 10 ksT, either near the 
mean-field critical point or the spinodal curve, the dis- 
tinction between nucleation and spinodal decomposition 
will be lost and our approach will require modification. 

We are not aware of any experimental investiga- 
tions of the shape and size of nuclei in transitions 
between the lamellar and cylindrical phases in di- 
block copolymer melts. However, it is worth not- 
ing some experimental observations on related systems. 
Koizumi et al. used transmission electron microscopy 
to study a blend of poly(styrene-feZocfc-isoprene) and 
homopolystyrene, which macrophase-separated to form 
diblock-rich, lens-shaped droplets in a homopolymer-rich 
matrix Microphase separation into a hexagonal 

lattice of cylinders occurred in the droplets, with the 
cylinder axis aligned along the short axis of the lens-like 
macrodomain, which had an aspect ratio of about 3.5. 
Perpendicular to the cylinder orientation, the droplet had 
an approximately circular cross-section of about 20 cylin- 
ders in diameter. While their observations closely resem- 
ble our results for the critical droplet at Ja = 0.45, their 
system is quite different from ours. Firstly, they are look- 
ing at macrophase separation, rather than nucleation. 
Secondly, in their droplet the cylinders terminate at the 
interface with the disordered, homopolymer-rich phase, 
while ours merge continuously into a lamellar structure. 
Thus, while the relationship, if any, between these ex- 
periments and our theory remains to be understood, it 
is intriguing that droplets similar to those predicted here 
exist in nature. 

Balsara and coworkers have examined the evolution of 
cylindrical order in poly(styrene-6Zocfc-isoprene) follow- 
ing a quench from the disordered state |27L |28L l29l | . From 
their depolarized light-scattering data they inferred that 
their nuclei had an aspect ratio of about 4, and that nu- 
cleation was occurring with the cylinders oriented along 
the long axis of the droplet. While they observe an 
anisotropic droplet shape, their cylinders are oriented 
oppositely to those in the present theory, and those in 
the experiments of Koizumi et al. |48j. Of course, since 
the situations studied are all different, it is difficult to 
draw any conclusions by comparison. However, we have 
repeated the analysis described here for a droplet of cylin- 
der phase in a metastable disordered background and our 



preliminary results indicate that the cylinders still align 
along the short axis of the droplet. Again, the interfacial 
free-energy is lower when the droplet interface lies per- 
pendicular to the cylinder axis. This discrepancy with 
the conclusions of Refs. [23, needs to be better 

understood. 

Nonomura and Ohta recently performed two- 
dimensional simulations of the nucleation and growth of 
a droplet of lamellar phase from a metastable cylinder 
phase 12^ . Although this is the reverse of the transition 
considered here, comparison is possible since the present 
theory for the interfacial profile and free-energy does 
not refer to which phase is on the inside of the droplet. 
Their simulations, performed in what here is the x-y 
plane, observe the epitaxy we assume. Their droplet 
interface is relatively sharp, with a width on the order 
of 2-3 cylinder spacings. Their data suggest that the 
droplet has a hexagonal shape in the the x-y plane, 
although their larger droplets, on the order of 20 
cylinders across, appear more circular. In the absence 
of a precise definition of the droplet surface in the 
simulations, these observations about the droplet shape 
can only be considered qualitative. It appears that 
the critical droplet size in the simulations is less than 
1-2 cylinder spacings, since droplets of all sizes grow. 
This small critical droplet size, less than that found 
here, may occur because of different model parameters 
used in the simulations, or because the droplet was 
nucleated heterogeneously by creating a dislocation pair 
in the cylindrical order, thereby lowering the barrier to 
nucleation. 

We have demonstrated the crucial connection between 
anisotropy in the interfacial free-energy and droplet 
shape, and Fig. [51 demonstrates the importance of using 
the proper critical droplet shape to accurately calculate 
the nucleation barrier. Thus it is important to consider 
the robustness of our results for the critical droplet shape 
and nucleation barrier. We consider the class of models 
for diblock copolymers where the free-energy is written 
in terms of a Landau expansion in the order-parameter, 
and whose gradient structure occurs at quadratic order 
in the order-parameter. This class includes the Landau- 
Brazovskii model, the Leibler model with local cubic and 
quartic coefficients j^], and the Ohta-Kawasaki model 
[33j . If the single-mode and slowly- varying amplitude ap- 
proximations are applied to this class of models one can 
show that the resulting gradient structure will reduce to 
that of Eq. (|2.25|) . Since the critical droplet shape de- 
pends crucially on the gradient structure of the theory, 
this suggests, for weak segregation, that the lens-like crit- 
ical droplet shape calculated here is robust. At higher 
segregation the droplet shape may be modified from that 
calculated here. 

The need to use a variational ansatz to obtain the free- 
energy of a planar interface can be eliminated by solving 
the amplitude model, Eqs. H2.26|l and H2.27|l . numerically 
in three dimensions for a droplet of cylinder phase in a 
background of metastable lamellar phase. Near the phase 
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boundary, we expect no qualitative change in either the 
droplet shape or the behaviour of the nucleation barrier 
from that reported here. The overall magnitude of the 
nucleation barrier will, however, be somewhat reduced 
since the elimination of the variational ansatz will lead to 
a reduction in the calculated interfacial free-energy. This 
approach will also enable us to go beyond the constraints 
imposed by the Wulff construction and the classical nu- 
cleation theory. It will allow an investigation of the en- 
tire metastable lamellar region from the phase boundary 
to the spinodal. The nucleation barrier should approach 
zero and the critical droplet dimensions should diverge at 
the spinodal 0- As the spinodal is approached and the 
droplet interface becomes more diffuse the connection be- 
tween interfacial free-energy and droplet shape becomes 
less clear, thus the droplet shape may be modified near 
the spinodal. It may be possible to make contact with 
other theories that address the spinodal limit of order- 
order transitions USE 111^2, 13, 16]. 

In addition to exploring the limits of the theory 
through the extensions discussed above, we plan to apply 
the formalism to nucleation at the sphere /cylinder tran- 
sition. The role of anisotropic fluctuations during this 
transition has been studied experimentally in Ref. |49|. 
This transition has been studied theoretically in Refs. 

assuming uniform intermediate states, but the 
nucleation of compact droplets has not been examined. 
Finally, we would like to develop a deeper understand- 
ing of the orientation dependence of the interfacial free- 
energy in terms of the structure of the lamellar/cylinder 
interface and the conformations of block chains at these 
interfaces. 



VI. CONCLUSION 

We have examined the nucleation of a droplet of sta- 
ble cylinder phase from a metastable lamellar phase using 
the single-mode approximation to the Brazovskii model 
for diblock copolymer microphases. By employing a vari- 
ational ansatz for the droplet interfacial profile, we find 
an analytic expression for the interfacial free-energy of 
an interface of arbitrary orientation between cylinders 
and lamellae. The interfacial free-energy is anisotropic, 
and is lower when the cylinder axis is perpendicular to 
the interface than when the cylinders lie along the inter- 
face. Furthermore, the interfacial free-energy is slightly 
lower when the lamellae are parallel to the interface 
compared to the perpendicular alignment. The droplet 
shape computed via the Wulff construction is lens-like, 
being flattened along the axis of the cylinders. As the 
mean-field critical point is approached along the lamel- 
lar/cylinder phase boundary the droplet becomes more 



flattened along the cylinder axis. We apply this informa- 
tion to compute the size of the critical droplet and the 
nucleation barrier within classical nucleation theory. We 
are able to make specific predictions for these quantities 
by connecting the phenomenological Brazovskii model to 
the many-chain Edwards Hamiltonian for diblock copoly- 
mers. The general trend is for the nucleation barrier 
to decrease and the critical droplet size to increase as 
the mean-field critical point is approached. The nucle- 
ation barrier is significantly reduced when the critical 
droplet shape is anisotropic instead of spherical, indicat- 
ing the importance of using the proper droplet shape. 
The theory should be valid near the lamellar/cylinder 
phase boundary from the lamellar/gyroid/cylinder triple 
point at /a ~ 0.45 to near the mean-field critical point. 
In this regime, droplets of size 30-400 cylinders across 
with aspect ratios of 4-10 and nucleation barriers of 30- 
40 ksT are typically found. Close to the mean-field criti- 
cal point fluctuations may modify this mean-field picture. 
Due to the variational approximation, the computed in- 
terfacial free-energies are upper bounds on the true free- 
energies, implying that the computed critical droplet di- 
mensions are also upper bounds on the true dimensions. 

Our theory makes specific predictions about the nucle- 
ation barrier and the size and shape of the critical droplet 
for nucleation of the cylindrical phase from metastable 
lamellae in diblock copolymer melts. We systematically 
calculate the interfacial free-energy for interfaces of ar- 
bitrary orientation between lamellae and cylinders. This 
work is an important first step towards a more sophis- 
ticated theory of nucleation in this system. The size of 
the critical droplets we find already suggests that a di- 
rect numerical attack on the nucleation problem, using 
Eq. (|2.1() . will be challenging. Although our focus has 
been on nucleation in diblock copolymers, this approach 
should work for any system in the Landau-Brazovskii 
class, given the appropriate model parameters. The nu- 
cleation scenario described here should be observable ex- 
perimentally, however experiments which study droplet 
nucleation along the lamellar/cy Under phase boundary 
in diblock copolymer melts have yet to be performed. 
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